Conditions for magnetically induced singlet d-wave superconductivity on the square 

lattice 



O 

o 

> 
O 



3 



o 

> 
On 
CO 

o 



o 
o 



X 



S. R. Hassan \ B. Davoudi 1 * 2 , B. Kyung 1 and A.-M.S. Tremblay 1 
1 Department de Physique and RQMP, Universite de Sherbrooke, Sherbrooke, Quebec, Canada J IK 2R1 
2 Institute for Studies in Theoretical Physics and Mathematics, Tehran 19395-5531, Iran 

It is expected that at weak to intermediate coupling, d-wave superconductivity can be induced 
by antiferromagnetic fluctuations. However, one needs to clarify the role of Fermi surface topology, 
density of states, pseudogap, and wave vector of the magnetic fluctuations on the nature and strength 
of the induced d-wave state. To this end, we study the generalized phase diagram of the two- 
dimensional half- filled Hubbard model as a function of interaction strength U/t, frustration induced 
by second-order hopping t'/t, and temperature T/t. In experiment, U/t and t! jt can be controlled by 
pressure. We use the two-particle self-consistent approach (TPSC), valid from weak to intermediate 
coupling. We first calculate as a function of t'/t and U/t the temperature and wave vector at which 
the spin response function begins to grow exponentially. D-wave superconductivity in a half-filled 
band can be induced by such magnetic fluctuations at weak to intermediate coupling, but only if 
they are near commensurate wave vectors and not too close to perfect nesting conditions where the 
pseudogap becomes detrimental to superconductivity. For given U/t there is thus an optimal value 
of frustration t'/t where the superconducting T c is maximum. The non-interacting density of states 
plays little role. The symmetry d x 2_ y 2 vs d xy of the superconducting order parameter depends on 
the wave vector of the underlying magnetic fluctuations in a way that can be understood qualitatively 
from simple arguments. 

PACS numbers: 74.20.Mn, 71.10.Fd, 05.30.Fk, 74.20.Rp 



I. INTRODUCTION 

Even before the discovery of high-temperature super- 
conductivity, research in organic and heavy fermion com- 
pounds lead to the suggestion that antiferromagnetic 
spin fluctuations can promote superconductivity with an 
order parameter that changes sign along the diagonal 
of the Brillouin zone, so-called d-wave superconductiv- 
ityL2£&. This question can be studied using the Hubbard 
model. At strong coupling, namely interaction strength 
U larger than the bandwidth, mean-field factorization of 
the equivalent t — J model already reveals the possibility 
of d-wave pairing. The situation is quite different at weak 
coupling. There, no mean-field factorization of the Hub- 
bard model leads to a d-wave superconducting state. In- 
stead, magnetic fluctuations are dominant. Pairing must 
then be seen in a two-stage process. First, sufficiently 
strong magnetic fluctuations are formed, then the corre- 
sponding low energy bosons can act as a glue for pairs, 
in a manner analogous to what phonons do in the stan- 
dard BCS theory. In other words, at strong coupling it 
seems that superexchange suffices to form pairs, while 
at weak coupling there is an intermediate process where 
mediating bosons must be formed before they can bind 
pairs 5 . 

There is now ample numerical and analytical evidence 
for d-wave superconductivity in the one-band Hubbard 
mo del&>LS^lLi2^i^S^Ll&^^ even though there 
are some dissenting voices 2 ^. While the coupling in high- 
temperature superconductors or the layered BEDT or- 
ganics is strong enough to lead to a Mott insulating state 
at half-filling, it is still interesting to understand more 
deeply the physics at weak to intermediate coupling in 



the absence of the Mott insulating state. There, numer- 
ical methods are limited by the fact that small system 
sizes do not allow magnetic fluctuations at incommensu- 
rate wave vectors. 

In this paper, we study the two-dimensional Hubbard 
model at half-filling. In two-dimensions, the Mermin- 
Wagner theorem implies that there is no long-range mag- 
netic order at finite temperature. Nevertheless, there is a 
crossover temperature at which the magnetic correlation 
length starts to grow exponentially as temperature de- 
creases. One enters the so-called renormalized classical 
regime. We thus study the influence of second-neighbor 
hopping t' and of interaction strength U on the magnetic 
crossover temperature. Berezinsky-Kosterlitz-Thouless 
superconductivity is allowed in two-dimensions and can 
be induced by these magnetic fluctuations. We answer 
the question of the relation between the wave vector of 
the dominant magnetic fluctuations and the symmetry 
and magnitude of d-wave superconductivity^ 2 -, as well 
as the question of the relative importance of the non- 
interacting density of states, the double occupancy, the 
magnetic correlation length, and the detailed shape of 
the Fermi surface on superconductivity. 

We will use the Two-particle self-consistent (TPSC) 
approach 2 ^!^^ that is non-perturbative and valid 
at weak to intermediate coupling. This method has 
been benchmarked by comparisons with Quantum Monte 
Carlo calculations wherever possible. It gives results for 
the pseudogap that are in agreement with experiments in 
electron-underdoped cuprate a 27 ' 28 . In addition, the su- 
perconducting transition temperature obtained 2 ^ by this 
method is very close to more accurate results from the 
Dynamical Cluster Approximation 12 . 
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We first introduce the model and the TPSC method, 
explaining in detail why it is a method of choice for this 
problem. We next present our numerical results, discuss 
various physical aspects that can be found in the t — t' — U 
model and compare our results with other results in ear- 
lier literature. We then discuss the relationship between 
d-wave superconductivity and magnetic fluctuations in 
terms of a simple physical picture and conclude. 

II. MODEL AND FORMALISM 



and 

< n 2 > + < n\ > -2 < n ini >= — ^ X s P (q), (3) 

where /3 = i,n =< rif > + < rij >,q = (q,iq n ) with 
q the wave vectors of an N-site lattice, and with iq n — 
2irinT the bosonic Matsubara frequencies. The Pauli 
principle < rtr a >=< n a > applied to the left-hand side 
of both equations with the TPSC expressions for Xsp, Xch 
on the right-hand side leads to 



We study the Hubbard Hamiltonian, 

(id),* ((id)),* 



(1) 



where Ci CT (q ) are annihilation (creation) operator for 
electrons of spin a at site i. The hopping matrix ele- 
ment t is between nearest-neighbors < ij > whereas t 
is between next- nearest- neighbors << ij >>, n- ltT is the 
density operator at site i, and U is the on-site interac- 
tion. Particle-hole transformation with and without a 
phase factor shows that at half-filling the results do not 
depend on the sign of either t or t' . 

The TPSC approach that we use in this paper satis- 
fies the Pauli principle < ro 2 >=< n a >, conservation of 
spin and charge as well as the Mermin- Wagner theorem. 
In addition, it includes renormalization of interaction 
(Kanamori-Bruckner) and has a self-energy that includes 
the vertex corrections that must be taken into account 
in the absence of the conditions necessary to satisfy Mid- 
gal's theorem. TPSC also achieves consistency between 
one- and two-particle quantities in the form of the sum 
rule relating potential energy to the trace of self-energy 
times Green function. At weak to intermediate coupling 
(U < 6t), detailed comparisons with benchmark Quan- 
tum Monte Carlo simulation o 13 ' 24 i 26 reveal that TPSC is 
accurate to a few percent. 

In this approach, spin and charge susceptiblities Xsp, 
Xch have random phase approximation (RPA)-likc forms 
but with two different effective interactions U sp and U c h 
that are then determined self-consistently. Although the 
susceptiblities have a RPA functional form, the physical 
properties of the theory are very different from the RPA 
because of the self-consistency condition on U sp and U c h- 
The necessity to have two different effective interactions 
for spin and charge is dictated by the Pauli exclusion 
principle < n 2 >=< n a > which implies that the lo- 
cal values of Xsp and Xch are related to only one local 
pair correlation function < riffii >. Indeed, using the 
fluctuation-disspation theorem in Matsubara formalism, 
we have the exact sum rules 

< n 2 > + < n\ > +2 < n T n x > -n 2 = — ^ Xch(q) 
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(5) 



with xo (q) the susceptiblities calculated from Green func- 
tions whose constant self-energy can be absorbed in 
the chemical potential. Hence, xo(q) takes the non- 
interacting form. 

If < > is known, U sp and U c h are determined 

from the above equations. The key quantity < n^ni > 
can be obtained from Monto Carlo simulations or by 
other means. However, it can also be obtained self- 
consistently by adding to above set of equations the re- 
lation 



U sp =g n (0)U, 5n (0) = 



< rif >< ni > 



(6) 



Equations (5) and (6) define a set of self-consistent 
equations for U sp that involve only two-particle quan- 
tities. U c h is determined by substituting in Eq.(4) the 
value of < n^ni > that is obtained from Eqs. (5) and 
(6). 

Following the procedure introduced in Ref. 25.29| one 
can show that the self-energy can be calculated from 
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+U chX c h (q)}G {1) Ck + q) 



(7) 



where the Matsubara frequency associated with k is 
fermionic and that associated with q is bosonic. In G^ 1 ' 
the self-energy is a constant that can be absorbed in the 
chemical potential so that effectively is the non- 

interacting Green function for half-filling. 

As usual we can compute the spectral function from 
A(k,u) = -ImG (2) (k, u)/ir where 



G ( V R (k,u) 



1 



iS — ek + /i 



E*(k j£ 



(8) 



The interacting chemical potential fi is found by fixing 
the number of particles calculated from G^ 2 ' R to half- 
filling. 
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The above formalism can be extended^ to compute 
pairing correlations. Basically, the above steps are re- 
peated in the presence of an infinitesimal external pair- 
ing field that is eventually set to zero at the end of the 
calculation. This allows us to obtain the particle-particle 
irreducible vertex in Nambu space from the functional 
derivative of the off-diagonal X^ 2 ) with respect to the 
off-diagonal Green function. The <i-wave susceptibility 
is defined by \d — Jjf dr (T t A(t) A^) with the <i-wave 



In this equation, <?<i(k) = (cos k x — cos k y ) for d x i_ y i 
symmetry and <7d(k) = sin k x sin k y for d xy symmetry. 
Details of the calculations! procedure are explained in 
Ref. H^. The only difference is that, in the present paper, 
double occupancy is always computed self-consistently 
with Eqs.(5) and (6). At high temperature the pair- 
ing susceptibility obtained from the above approach is in 
quantitative agreement with benchmark Quantum Monte 
Carlo simulations 2 ^. However, because of the sign prob- 
lem, the latter are not available at very low tempera- 
tures. Nevertheless, temperatures in the simulations are 
low enough to observe the effect of the pseudogap. 

Since the expression for the susceptibility Eq. (JS|) con- 
tains only the first two terms of the infinite Bethe- 
Salpeter series, we use the temperature where the con- 
tribution of the vertex part (exchange of one spin and 
charge fluctuation) becomes equal to that of the direct 
part of the <i-wave pairing susceptibility 2 ^ as a rough es- 
timate for the transition temperature for <i-wave super- 
conductivity. In other words, we look for the equality 
of the sign and the magnitude of the two terms appear- 
ing in Eq. ©. This choice is motivated by the state- 
ment that 1/(1 — x) ~ 1 + x diverges when x = 1. The 
true Kosterlitz-Thouless transition temperature in two 
dimensions is expected to occur below the temperature 
determined from the Bethe-Salpeter equation. The tem- 
peratures that we will find for the transition temperature 
are thus upper bounds. 

Strictly speaking, TPSC looses its quantitative accu- 
racy when it is taken below the crossover to the renor- 
malized classical regime. In cases where we will find the 
superconducting T c far below the crossover to the renor- 
malized classical regime, we expect our results to be only 
qualitatively correct. Nevertheless, as a further check 
of the validity of this approach, we point out that de- 
tailed finite size studies near optimal doping (x = 0.10) 



pair creation operator At = J2iJ2 7 9 il) 4-[4+ 7 V 
being a form factor for a representative gap symmetry, 
(3 = 1/T, T T is the time-ordering operator, and r is imag- 
inary time. The final expression for the d-wave suscep- 
tibility in the zero-frequency limit can be written as fol- 
lows, if we use the notation where k and k' denote both 
wave vector and fcrmionic Matsubara frequency, 



(9) 



using the dynamical cluster approximation^ 2 - find that 
T c = 0.023£ when U = 4t, t' = while with TPSCS the 
corresponding value is T c = 0.029i. 



III. NUMERICAL RESULTS AND 
COMPARISONS WITH OTHER APPROACHES 

Let us first consider magnetic properties. In two di- 
mensions, there is no long-range magnetic order at fi- 
nite temperature because of the Mermin- Wagner theo- 
rem. Nevertheless, there can be a crossover temperature 
Tx where the antiferromagnetic correlation length (£) be- 
gins to increase exponentially. In such a case, £ becomes 
infinite only at zero temperature where long-range order 
is established. One can define the crossover tempera- 
ture Tx as the temperature at which the antiferromag- 
netic correlation length £ becomes equal to the thermal 
de Broglie wavelength. At this temperature, a pseudo- 
gap begins to appear at hot spots in the zero frequency 
single-particle spectral weight^. Hot spots are points 
of the Fermi surface that are linked by the wave vector 
where the magnetic fluctuations are becoming large. 

In practice, as soon as the correlation length becomes 
relatively large compared with the lattice spacing, the 
rise in £ as temperature is lowered is rather sharp. For 
convenience, We thus choose Tx as that temperature 
where Xspilx, q y , 0)/xo(q x , q y , 0) = 200, with q x and q y 
the wave vectors where the response function has a max- 
imum. We chose 200 for the sake of computational time 
and checked that the result did not change appreciably 
if we chose a higher value. 

In the top panel of Fig. [1] we show the crossover dia- 
gram of the frustrated Hubbard model on the square lat- 
tice in the U—t' plane without superconductivity. All pa- 
rameters with dimension of energy are measured in units 
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FIG. 1: (Color online) Crossover diagram at n = 1 in the U—t 1 
plane, where U and t' are measured in units of hopping t. In 
the top panel there is no superconductivity. In the bottom 
panel, the boundaries where d-wave superconductivity can 
appear are indicated. The color scale represents the value 
of the crossover temperature for the spin response function. 
ISDW refers to incommensurate spin-density wave and PM 
to paramagnetic phase. 



of nearest-neighbor hopping t. Boltzmann's constant is 
also taken as unity. The color code in Fig.l represents 
the variation of the crossover temperature. The areas 
indicated with (tt, tt), ISDW, (0, tt — 6), and PM are the 
commensurate, incommensurate, incommensurate close 
to (0, tt) (S approaches as t' increases in that region (see 
Fig. 2)) and nonmagnetic (paramagnetic) metallic state, 
respectively. The PM state indicates that the peak of the 
spin response function does not grow even at temperature 
as low as T = 0.01. 



We note that at small t', the largest peak value of 
the spin response function is at the (tt, tt) wave vector, 
which is expected to mediate d x i_ y i superconductivity. 
For t' > 0.8 on the other hand, the response is peaked at 
the (0, tt — S) wave vector, expected to be relevant for d xy 
superconductivity. Due to the exponentially increasing 
correlation length with decreasing temperature, a small 
coupling between two dimensional planes would lead to 
long range spin order at finite temperature, especially for 
small and t' > 0.8. 

Our crossover diagram (without superconductivity) 
can be compared with zero-temperature results obtained 
by other methods. For t' < 0.5 there are results of 
the optimized variational Monte Carlo (VMC) method 
in Ref. |30j. For U < 6, their antiferromagnetic region 
extends only up to \t'/t\ < 0.2, a result quite different 
from other studies including ours. For t' < 1 the phase 
diagram in the U—t' plane has been obtained by the path- 
integeral renormalization group (PIRG) approach 31 . Re- 
stricting ourselves to the region U < 6 of that paper, we 
find that the area of our PM state is smaller. However, 
the region where (tt, tt) fluctuations occur is quite close. 
They find a (tt, 0) phase for t' > 0.7 like us but at U > 8 
Due to the finite size of the PIRG calculations, they are 
not sensitive to the incommensurate spin density waves 
that we find. That may explain why their paramagnetic 
region is larger. Also, in the region where we find ISDW 
they find, at large enough U, a non-magnetic Mott insu- 
lator (NMI) and at even larger U a (tt/2,tt) phase. Our 
TPSC approach cannot describe the Mott transition, so 
the question is whether that transition can occur at val- 
ues as small as about U = 4 as found in Ref. |3jJ. That 
seems to be possible within the Variational cluster ap- 
proximation (VCA)22, but VCA overestimates the effect 
of interactions. The (tt, tt) region in the latter work is 
larger than ours, again perhaps because the small sys- 
tem sizes are not sensitive to incommensurate density 
waves. However, the (0, tt — S) region is quite close to 
ours. We should also keep in mind that we are quoting 
the position of the maximum spin susceptibility at the 
crossover temperature. That position may change when 
zero temperature is reached. In fact, quite generally, it 
is expected that the magnetic fluctuations become more 
incommensurate as temperature decreases. 

We open a parenthesis on the incommensurability dis- 
cussed above. It can be contrasted with what is found in 
spin models that correspond to the large U limit of the 
Hubbard model. Fig.[5]shows the classical ground state 33 
of the Ji(= t 2 /4U)-J 2 (= (t') 2 /4U) Heisenberg model on 
the square lattice as a function of t'/t — yj J\j 'J 2 . Com- 
mensurate (tt, tt) ordering is found for t' < *J\t while 

for t' > \J\t two independent AF sublattices appear as 
the ground state. Thermal or quantum fluctuations se- 
lect a collinear phased with (tt, 0) or (0, tt), and (q x ,^) 

and (TT,q y ) for t' — We can see in Figs. [3] and [4] 

that at weak to intermediate coupling, Fermi surface ef- 
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FIG. 2: The classical ground state- 3 of Ji — J2 Heisenberg 
model at T=0 as a function of t' obtained by minimizing the 
dispersion relation J(q) = 2[cos(q x ) + cos(q y )\ + 2t /2 [cos(q x + 
q y ) + cos(q x - q y )]. 
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^(0.5 - k™ ax ) 2 + (0.5 - k™ ax ) 2 
typical position of the wave vector where the maximum spin 
response function is located at Tx as a function of t' for 
U — 2, 4 and 6 all at n — 1. All energies are in units of t. 



fects as well as thermal and quantum fluctuations reduce 
the region where (tt, tt) commensurate ordering is found, 
relative to the classical strong coupling theory (U >> t). 

We now turn to the regions, illustrated in the lower 
part of Fig.l, where spin fluctuations can induce d-wave 
superconductivity. There are only five real irreducible 
representations of the tetragonal group that correspond 
to singlet pairing 3 ^. Amongst these, only B lg (d x 2_ y i) 
and £?2g (d xy ) have an internal structure that can take ad- 
vantage of spin fluctuations in the two-dimensional plane 
(see discussion section below). Incommensurate fluctua- 
tions could conceivably lead to triplet superconductivity 
in one of the five odd irreducible representations, a pos- 
sibility we do not consider here. 

At n = 1 and t' = 0, it was shown in Ref. that 
strong antiferromagnetic fluctuations create a pseudo- 
gap over the perfectly nested Fermi surface, suppressing 
the possibility of d-wave superconductivity. By including 
frustration in the form of second-neighbor hopping, the 
pseudogap becomes limited to the hot spot regions^, al- 
lowing the d x 2_ y 2 superconducting state to appear. As 
the ordering wave vector becomes more incommensurate 
with increasing t', the conditions necessary for d x 2_ y 2 
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FIG. 4: (Color online) The position of the wave vector where 
the maximum spin fluctuations occur in the Brillouin zone for 
n = 1, U = 4 as a function of t'. Numbers near data points 
indicate the value of t' . 



superconductivity are no-longer satisfied (see discussion 
section below). This explains the finite region over which 
d x 2_ y 2 appears in Fig. [TJ Magnetic fluctuations near 
(tt, 0) lead to d xy symmetry (see discussion section). The 
area between the solid lines with circles (green) and the 
dashed lines with triangles (blue) indicate respectively 
the d x 2_ y 2 and d xy superconducting regions. We can see 
that the window of d-wave superconductivity decreases 
upon decreasing U and vanishes at (U = 1.9, t' — 0.09) 
for d x 2_ y 2 and at (U = 3,<' = 1.85) for the d xy state. 
When U is not large enough in the presence of frustra- 
tion, antiferromagnetic fluctuations are not large enough 
to lead to d-wave superconductivity. The large U(> 6) 
effect on the window where d-wave superconductivity ex- 
ists cannot be captured within TPSC since this approach 
is valid only from weak to moderate coupling. 

The region t' < 0.5 in the above phase diagram bears 
some similarity with that found by T = extrapolation of 
the fluctuation exchange approximation (FLEX) 36 . The 
value of (U, t') at which the d x 2_ y 2 superconducting state 
ceases to exist is around (3.38, 0.33) in that approach. 
The d x 2_ y 2 window increases rapidly with increasing U 
and becomes approximately constant beyond U = 4.0. 
(We are comparing the results of this paper for positive 
t' since they do not show the expected symmetry under 
change of sign of t' at n = 1.) In the T = VCA study 3 ^ 
mentioned above, the d x 2_ y 2 state appears all the way 
down to U — and up to t' — 1. Our results are ex- 
pected to be more accurate in the weak-coupling region. 
With TPSC we cannot address the question of T = 
homogeneous coexistence of superconductivity and anti- 
ferromagnctism (superconducting antiferromagnet). 

The d x 2_ y 2 superconducting critical temperature T c as 
a function of t' at U = 2.5, 3 and 4 is depicted in Fig. [5] 
The inset shows T c (t') for the d xy case at U = 3.6,4. 
The maximum T c increases with increasing U, by con- 
trast with what is found at strong couplin g 11 ^ 6 ' 37 . T c 
for the d x 2_ y 2 state as a function of t' at fixed U shows a 
behaviour that is similar to that found earlier as a func- 
tion of doping 26 for t' = 0. In fact, even the maximum 
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FIG. 5: (Color online) The d x 2_ y 2 superconducting critical 
temperature T c as a function of t' at U = 2.5, 3 and 4 for 
n — 1 . The inset shows the d xy superconducting critical 
temperature T c as a function of if for U = 3.6 and 4. 



FIG. 7: (Color online) g = gn(0) as a function of t' at U = 4 
and 6 for n = 1 at the crossover temperature Tx ■ Inset shows 
Tx as a function of t' at U = 4 and 6 for n = 1. 
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FIG. 6: (Color online) Logarithm base ten of the antifer- 
romagnetic correlation length (in units of the lattice spac- 
ing) as a function of inverse temperature for three values of 
t' = 0.15, 0.21, 0.31 at U = 4 for n = 1. The value of T c for 
the corresponding t' is shown on the plot. 



value of T c that can be reached for given U is compa- 
rable in the two cases. There are some similarities in 
the physics. The antiferromagnetic fluctuations need to 
be large enough and also at the correct wave vector to 
lead to superconductivity, but if they are too large with 
hot spots that cover the whole Fermi surface, supercon- 
ductivity disappears. In FLEX 36 there is no pseudogap 
but the T = extrapolation leads to antiferromagnetism 
destroying superconductivity at t' = 0. 

We now study in turn how incommensuration, anti- 
ferromagnetic correlation length, double occupancy and 
Fermi surface effects influence d-wave superconductivity. 

Fig. [3] illustrates at the crossover temperature Tx 
the incommensurability of the wave vector where the 



maximum spin response occurs. The quantity f(t') = 

^(0.5 - k™ ax ) 2 + (0.5 - k^ ax ) 2 is plotted as a function 

of t' at U = 2, 4 and 6. Here, k™ ax and k™ ax are the 
x and y components of the wave vector where the spin 
response function is maximum. In f(t'), k™ ax and k™ ax 
are measured in units of 2ir so that 0.5 corresponds to rr. 
The position of the peak is at (0.5, 0.5) for small t' and 
approaches (0.5,0) with increasing t'. The peak position 
at Tx does not strongly depend on U even though Tx de- 
pends on U. To be more explicit, we show in Fig. 3] how 
the wave vector moves in the Brillouin zone. Since the 
spin correlation function has typically several maxima, 
the highest peak determining the largest spin correlations 
sometimes suddenly jumps from one place to another in 
momentum space, which clearly happens near t' = 0.55 
and 0.7. Since the Fermi surface topology changes near 
if = 0.71, this is not surprising. 

As can be seen by comparing Fig. |4] and Fig. [5l the 
maximum value of T c for d x i_ y i superconductivity al- 
ways occurs in the region where the antiferromagnetic 
fluctuations are commensurate at (7r,7r), namely for < 
if < 0.32. Similarly, the maximum T c for d xy supercon- 
ductivity occurs where magnetic fluctuations are nearly 
commensurate at (0,7r — S), although in the latter case, 
strict commensurability is never achieved. The maximum 
T c at U = 4, for example, is smaller for d xy compared 
with d x 2_ y 2 . 

The antiferromagnetic correlation length £ is a mea- 
sure of how strong are the antiferromagnetic fluctua- 
tions. Here, we find £ by calculating the curvature of 
the zero-frequency interacting spin susceptibility at the 
maximum. To assess the importance of this factor on 
superconductivity, the logarithm of £ is shown in Fig. [6] 
for three different values of if at U = 4. An arrow in- 
dicates for what length we find a superconducting T c . 
The antiferromagnetic correlation length at T c is small- 
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FIG. 8: (Color online) Non-interacting Fermi surfaces (black 
solid lines) at half-filling and wave vectors of the dominant 
spin fluctuations for four values of if and < k x . y < 2n. For 
negative values of if the plots are identical if — n < k x ,y < 7r. 
The arrows represent the wave vectors that lead to the largest 
spin susceptibility, a) if = b) t' = 0.3. Dashed line is the 
antiferromagnetic Brillouin zone, c) t' = 1.2. d) if = 2.0. 
In the last two plots, the dashed lines are separated by (0, 7r) 
and (tt, 0). 



est for the largest T c . Near perfect nesting (smallest t'), 
where pseudogap effects are important, superconductiv- 
ity is hindered by the removal of states near zero en- 
ergy^. To compensate, the antiferromagnetic correlation 
length has to become very large before superconductiv- 
ity can set in, as we can see from t' = 0.15 on the figure. 
We can also compare at T c the thermal de Broglie wave- 
length vp/irT in units of the lattice spacing (and ks = 1, 
h = 1) with the value of £. In our units, the thermal 
length takes a value close to 1/T on the horizontal axis 
of Fig [HI It is only for the largest value if = 0.31 that 
the superconducting T c occurs well above the renormal- 
ized classical regime, i.e. for the thermal length much 
larger than £. The important qualitative fact is that T c 
occurs below the renormalized classical regime Tx when 
its occurrence is hindered by the pseudogap and that it 
can occur above Tx otherwise. Near the maximum T c , 
the two lengths are of the same order, in other words, for 
that t' = 0.21, T c and Tx are of the same order. But in 
all cases £ is much larger than the lattice spacing. 

For large t', where d xy superconductivity occurs, the 
antiferromagnetic correlation length is large, correspond- 
ing to the smaller T c obtained for that pairing symmetry 
(not shown). 

Returning to the factors influencing d-wave supercon- 
ductivity, it is meaningful to ask how double occupancy 
in the underlying normal state influences the appear- 
ance of superconductivity. One might naively think that 



larger double occupancy would be favorable to super- 
conductivity. A counter argument is that larger interac- 
tion strength, hence smaller double occupancy, is more 
favourable to d-wave superconductivity, hence should 
correspond to larger T c . We show in Fig. [7] double oc- 
cupancy normalized to its Hartree-Fock value, namely 
g = gn(0), as a function of t' at U = 4 and 6 at the 
crossover temperature Tx- The value of Tx as a func- 
tion of if for the same values of U is shown in the in- 
set. The maximum in double occupancy in the region 
0.4 < if < 0.8 corresponds to weaker spin fluctuations 
(as expected from the sum rule Eq.(5)) and vanishing 
T c . This is in agreement with what one can observe (yel- 
low region) in Fig. [TJ Comparing the overall values of 
g = 3ix(0)j we check, as expected, that double occupancy 
is smaller at larger U. The smaller double occupancy for 
larger U corresponds to a larger value of the maximum 
T c [if) in Fig. Hence, from all of the above one might 
conclude that smaller double occupancy is better for d- 
wave superconductivity. However, for the d x 2_ y 2 case 
the maximum in T c as a function of t' occurs for both 
values of U in a region where double occupancy is con- 
stant, independent of t' . Hence, one cannot conclude that 
smaller double occupancy is more favorable to supercon- 
ductivity. In fact the opposite trend is seen in the fact 
that the maximum T c at U = 4 in the d xy case occurs 
in a region where double-occupancy is large. Even more 
striking is the fact that the value of double occupancy in 
the d xy case (t' > 1) is monotonically increasing with t' 
while T c is not. Clearly, double occupancy, a local quan- 
tity, is not directly correlated with the value of T c for 
order parameters with nodes at the origin. 

We move to the influence of the shape of the Fermi 
surface on d-wave superconductivity. First of all, the 
shape of the Fermi surface influences the wave vector at 
which magnetic fluctuations are largest. Fig.[5]illustrates 
the fact that at t' = perfect nesting occurs for (%, tt) 
and symmetry related vectors while for if = 2 one has 
nearly perfect nesting at (0, 7t — 5). The change in Fermi 
surface topology occurs around if — 0.71. There is no 
further change in topology between if — 0.71 and if = 2. 
The symmetry of the induced superconductivity is tied 
to the wave vectors of the magnetic fluctuations, namely 
d x 2_ y 2 for (tt, it) and d xy for (0, 7r). No singlet supercon- 
ductivity occurs in the vicinity of the change in Fermi 
surface topology because of the incommensurability and 
smallness of the fluctuations. 

The single particle density of states of the non- 
interacting system is not a major factor in the occur- 
rence of d-wave superconductivity. This is demonstrated 
in Fig. [9] that shows that there is no correlation with the 
maximum T c in Fig. O and peaks in the single particle 
density of states. There is a jump in density of states 
near t' = 0.71 where Fermi surface topology changes. 
The only case where the Fermi surface coincides with 
a van Hove singularity in the single-particle density of 
states is at if = and in that case there is no d-wave su- 
perconductivity, contrary to naive expectations. This is 
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FIG. 9: (Color online) Non-interacting single particle density 
of states at the Fermi level No as a function of t' , calculated 
at half-filling with an energy resolution 0.02. Spin degeneracy 
is included. 




FIG. 10: (Color online) Illustration of the scattering by 
(— 7r,7r) fluctuations for d x 2_ y 2-~w&ve superconductivity and 
by (—tv + S, 0) for d xy -wa,ve superconductivity. 



because the nesting allows strong antiferromagnetic fluc- 
tuations at (tt, 7t) to open a pseudogap over the whole 
Fermi surface. This is illustrated in Ref. HUM The 
pseudogap removes density of states from the Fermi en- 
ergy and is detrimental to d-wave superconductivity 26 . 
The same phenomenon is seen at t' = 2 where there is 
also nearly perfect nesting. In the presence of a finite t' , 
the pseudogap is created only at hot spots and d-wave 
superconductivity can arise. 



IV. DISCUSSION 

Our estimate of T c from the normal state rests on the 
fact that at T c the vertex equation deduced from the 
Bethe-Salpeter equation becomes singular. In that case, 
the vertex equation has the same structure, symmetry- 
wise, as the BCS gap equation. One can thus understand 
very simply the relation between the symmetry of the d- 
wave order parameter and the underlying magnetic fluc- 
tuations. Normally, one would project the interaction 



potential in the different pairing channels^^. A less 
rigorous but more intuitive approach follows Ref. and 
starts from the BCS gap equation 

A P = -i/^V(p-p')^(l-2n(£; p ,)). (10) 

where E p / is the quasiparticle energy and n(E p i) the 
Fermi function. Even though the effective potential 
V (p — p') for singlet pairing due to spin fluctuations 
is always positive (repulsive), it can lead to an effective 
negative (attractive) potential for d 2 .2_ y 2-wave supercon- 
ductivity if it is peaked at (— 7T, it) and symmetry related 
wave vectors, as illustrated on the left of Fig. [TO] In- 
deed, take a value p along the x axis where the gap is 
positive. Then the largest contribution to the integral 
over p' will come from a region around p' located at 
±90°, where the gap is negative, such that the condi- 
tion p p' = (—tt, ±7r) is satisfied as closely as possi- 
ble. The change in sign of the gap allows a solution of 
the gap equation even with V (p — p') positive. (Clearly 
also, increasing the strength of the interaction will lead 
to a higher T c .) The analogous argument explains why 
(±7r, 0), (0, ±7r) scattering favors d xy -wsve superconduc- 
tivity, and why magnetic fluctuations that are peaked at 
incommensurate wave vectors are not favorable to any 
of the allowed singlet pairing symmetries on the square 
lattice. (In addition for a given U the magnitude of fluc- 
tuations at incommensurate vectors is smaller than at 
commensurate vectors). While we have not studied the 
following question in detail, we note that incommensu- 
rate wave vectors may be favorable to pairing in triplet 
channels such as p-wave, (E u )2^. There is however a 
change in sign of the effective interaction in the triplet 
channelii^. That kind of mechanism may occur in the 
heavy fermion compound CeCu2Si?r& for example. 

According to Anderson^, that physics whereby the 
pair state is orthogonal to the repulsive core of the 
Coulomb interaction was explained by Pitaevskii and 
Bruckne r 39 ' 40 . Note also that in our approach, we neglect 
the feedback of superconductivity on spin fluctuations 4 ^. 

As discussed in the previous section, in the present case 
the non-interacting density of states has little influence 
on the value of T c . Self-energy effects can strongly mod- 
ify the effective density of states by creating a pseudogap 
that is detrimental to superconductivity. Self-energy ef- 
fects are not so important in calculations that do not 
contain the pseudogap^. We saw the suppression of su- 
perconductivity by the pseudogap at t' = for example. 
That suppression is already apparent at high tempera- 
ture in Quantum Monte Carlo calculations^ 6 - and is repro- 
duced by TPSC within a few percent in that regime. In 
apparent contradiction with the suppression of supercon- 
ductivity by strong antiferromagnetic fluctuations in the 
presence of nesting, in quantum cluster approache s 12 ' 16 , 
one does not see the suppression of d-wave superconduc- 
tivity at half-filling for t' = at weak to intermediate cou- 
pling (i.e. below the Mott transition), unless long-range 
antiferromagnetism is allowed. This is because for the 



9 



small cluster sizes considered, the normal-state magnetic 
correlation length cannot become large enough to cre- 
ate a weak coupling pseudogap. Indeed, a weak coupling 
pseudogap can barely be seen in Ref. |43| for cluster sizes 
4x4, even though they are larger than the 2x2 clusters 
in Ref. O It is only at even larger cluster sizes that the 
weak coupling pseudogap appears in DC A calculations 43 . 
The difference between weak and strong coupling pseu- 
dogaps has been discussed for example in Refs. l4#45l . 

In electron-phonon mediated superconductivity, it has 
been found that the phonon frequencies near 7kgT are 
the most efficient ones for pairing 4 ^. In the present case 
the wave vector of the magnetic fluctuations is also im- 
portant. Although we have not studied the importance of 
the frequency dependence in detail, it is possible to make 
the following remarks. Take d x 2_ y 2 superconductivity. 
For small t', where superconductivity is reduced by self- 
energy effects coming from the pseudogap, the important 
frequencies are less than temperature in energy units be- 
cause T c occurs in the renormalized classical regime. At 
larger t', where T c is reduced by incommensuration ef- 
fects and deviations to nesting, the pseudogap occurs at 
lower temperature than T c hence finite frequencies larger 
than temperature are still important. Optimal T c occurs 
between these two cases. 



V. CONCLUSION 

We have studied the conditions for the appearance of 
magnetic-fluctuation induced d-wave superconductivity 
in the half-filled Hubbard model in two dimensions. We 
have shown that at weak to intermediate coupling, the 
symmetry of the d-wave order parameter is determined 
by the wave vector of the magnetic fluctuations. Those 
that are near (tt,tt) lead to d x 2_ y 2-wave (Bi g ) supercon- 
ductivity while those that are near (0, tt) induce d xy -wave 
(Big) superconductivity. The dominant wave vector for 
magnetic fluctuations is determined by the shape of the 
Fermi surface so we find that d x 2 _ y 2 -wave superconduc- 
tivity occurs for values of if that are relatively small while 
d^y-wave superconductivity occurs for t' > 1. For inter- 
mediate values of t' the magnetic fluctuations are smaller 



and incommensurate so no singlet superconductivity ap- 
pears. The maximum value that T c can take as a function 
of t' increases with interaction strength. All of the above 
can easily be understood physically from simple BCS- 
like arguments. However, contrary to what is expected 
from BCS, the non-interacting single-particle density of 
states does not play a dominant role. With Fermi sur- 
faces where hot spots can create a pseudogap on nearly 
all the Fermi surface (as in the t' = case), self-energy 
effects hinder d-wave superconductivity, even though the 
strength of magnetic fluctuations can be very large. In 
that case, the large inelastic scattering rates are pair 
breaking and remove states near the Fermi level, decreas- 
ing the tendency to d-wave superconductivity. There is 
thus an optimal value of t' (frustration) for superconduc- 
tivity. For d x 2_ y 2 superconductivity in underfrustrated 
systems (small if) T c occurs below the temperature Tx 
where the crossover to the renormalized classical regime 
occurs. In other words, at T c the antiferromagnetic corre- 
lation length is much larger than the thermal de Broglie 
wave length. The opposite relationship between these 
lengths occurs for overfrustrated systems (t' larger than 
optimal) where T c is larger than Tx- The two temper- 
atures are comparable for optimally frustrated systems. 
But in all cases, at T c the antiferromagnetic correlation 
length is larger than the lattice spacing. 

Comparisons with experiments on heavy-fermion com- 
pounds, quasi one-dimensional organic metals and the 
question of the interplay of incommensurate magnetic 
fluctuations and triplet pairing is left for future work. 
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